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Understanding the organization of collective motion in biological systems is an ongoing challenge. 
In this Paper we consider a minimal model of self-propelled particles with variable speed. Inspired 
by experimental data from schooling fish, we introduce a power-law dependency of the speed of each 
particle on the degree of polarization order in its neighborhood. We derive analytically a coarse- 
grained continuous approximation for this model and find that, while the variable speed rule does 
not change the details of the ordering transition leading to collective motion, it induces an inverse 
power-law correlation between the speed or the local polarization order and the local density. Using 
numerical simulations, we verify the range of validity of this continuous description and explore 
regimes beyond it. We discover, in disordered states close to the transition, a phase-segregated 
regime where most particles cluster into almost static groups surrounded by isolated high-speed 
particles. We argue that the mechanism responsible for this regime could be present in a wide range 
of collective motion dynamics. 



I. INTRODUCTION 

Bacterial colonies fT] , insect swarms P| , bird flocks [3] , 
and fish schools 4] are all examples of biological sys- 
tems that display distinct collective motion. While they 
differ in many specific aspects, they also share an impor- 
tant characteristic: Each one is made up of individual 
organisms that are self-propelling and move guided by 
interactions with their surrounding neighbors. When we 
strip down collectively moving systems in this way, we 
obviously leave out important biological details, but in 
return we obtain a good starting point for a theoretical 
understanding of collective motion. Our central concern 
in this paper is to introduce and study a very simple 
model (a minimal model) where individuals move with 
variable speed, and to show through theory and simula- 
tions that it displays new properties not found in models 
with constant speed. 

In recent years, research on self-propelled (or active) 
systems has grown steadily. Starting with the seminal 
work by Vicsek et al. ^5^ , we now find a range of theoret- 
ical, numerical and experimental studies [51115) concern- 
ing a number of aspects of their dynamics. In Vicsek's 
original work 0, a minimal model of self-propelled col- 
lective motion — now known as the Vicsek model — was in- 
troduced, similar in flavor to the equilibrium XY-model 
[16j . but exhibiting new physical properties due to its 
non-equilibrium nature. 

In the standard Vicsek model, each individual is de- 
scribed as a point particle with a given position and 
heading direction. All particles advance with the same 
constant speed and decide their next heading direction 
based on the headings of all particles in their local neigh- 
borhood. Particles are also subject to noise, which alters 
their chosen direction. Under these simple conditions, 
the Vicsek model exhibits a dynamical phase transition 
from a disordered state (with no polarization order) to 



an ordered state (where particles align and advance to- 
gether) as the noise level is decreased, or the density of 
particles increased. 

Following its initial formulation, many variations of 
the original Vicsek model [lOl [17] have been introduced 
to carry out analytical [51 [Tl] [T^ [TH] and numerical [51 [TU] 
studies of collective motion. However, almost all of these 
models consider particles that advance with equal con- 
stant speed. But in biological systems the speed of each 
individual could vary in response to its neighbors' dy- 
namics, and it is reasonable to expect this variation to 
be important for the resulting collective motion. Fur- 
thermore, in recent experiments on fish schooling, it has 
been observed that fish swimming in disordered regions 
typically move slower than those in ordered ones [5DJ [H] . 
This observation supports the idea that variable speed 
may play an important role in the collective dynamics of 
real groups of self-propelled individuals. 

In this Paper, we introduce an experimentally moti- 
vated minimal model describing the collective motion in 
groups of self-propelled particles with variable speed, and 
study it analytically and numerically. The model is al- 
most identical to the standard Vicsek model, but here 
both the speed and orientation of each particle depends 
on the state of its local neighborhood. While the choice 
of speed dynamics is motivated by our experimental re- 
sults on schooling fish, we emphasize that the model is 
not intended to function as a detailed replicate of this 
biological system. Note that models similar to the one 
presented here were used in two recent numerical studies 
focusing on how variable speed can enhance convergence 
to an ordered state [22]. These studies, however, did not 
address the analytical results or novel spatial dynamics 
discussed here. Starting from our model, we derive the 
corresponding hydrodynamic equations of motion for the 
density of particles and for the order parameter field, 
which undergoes a symmetry breaking transition at the 
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critical point, like in the original Vicsek model. Using 
these equations, we find an analytical relation between 
the coarse-grained particle speed and the local density 
of particles, which we verify numerically. We then study 
our variable speed model numerically, characterizing its 
ordered and disordered phases as a function of the noise, 
mean particle density, and variable speed rule. In partic- 
ular, we discover in the disordered regime a novel state 
where static cluster are formed, containing particles mov- 
ing at speeds close to zero. Within these clusters, each 
particle receives conflicting heading information from its 
neighbors and thus advances at low speeds, which in turn 
limits its ability to spread its own heading information 
throughout the group. We hypothesize that this mecha- 
nism could also be present in more realistic systems. 

The Paper is organized as follows. In Section II we in- 
troduce our variable speed model using new experimen- 
tal results to motivate the particle speed dynamics. In 
Section III we derive the corresponding coarse-grained 
hydrodynamic equations of motion and compare them to 
the constant speed case. Section IV presents numerical 
results that confirm the validity of the analytical approx- 
imations and explore regimes beyond these, where the 
hydrodynamic description is no longer valid. We study 
here the static clusters described above. Finally, Sections 
V and VI contain our discussion and conclusions. 



II. VARIABLE SPEED MODEL AND 
EXPERIMENTAL MOTIVATION 

A. Experimental background 

Fish schools are a clear example of systems exhibit- 
ing collective motion. To motivate the choice of speed 
dynamics in our variable speed model we examined ex- 
perimentally the relationship between individual speed 
and local polarization order in a school of 300 golden 
shiners [Notemigonus crysoleucas) swimming freely in a 
shallow tank. A snapshot from the experiment is shown 
in Fig.[lja). The state of each individual i at a given time 
t is defined by its position (t) (defined as the centroid of 
the fish's image), speed Vi{t) and heading direction unit 
vector fii (t) , determined from automated video tracking. 
Further details on the experimental setting, protocol and 
tracking system can be found in Katz et al. [Hj. 

To quantify the level of heading alignment in the local 
neighborhood of an individual fish i we first define the set 
Ui that contains all neighbors j within a given interaction 
radius f, such that | < r. Note that Ui also includes 

the focal fish i. We then define Xi a measure of the 
local polarization order around individual i, given by 
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jeUi 



(1) 



We analyzed the experimental data by finding for ev- 
ery individual fish in each of the 300 000 frames recorded 
its speed Vi{t) and local order Xi- We use here f « 15.5 
cm, but verified that similar results are obtained for dif- 
ferent values of f. This interaction zone is illustrated in 
Fig. [ija) by a white circle centered around a focal fish 
marked by a black dot. We then generated a histogram of 
the Vi{t), Xi data. For each bin value of x, we normalized 
the speed in order to get the probability distributions for 
the speed as a function of the local polarization order. 
The resulting data is visualized as a smoothed contour 
plot in Fig. [ijb). It is apparent that there is typically 
a strong (superlinear) relationship between the speed of 
an individual fish and the alignment of its local neigh- 
borhood. Note that this result does not reveal anything 
about the causality, since we only show here correlations 
within the data. These could be a consequence, for ex- 
ample, of having individuals that advance slower when 
they find themselves in a disordered environment or that 
are unable to align to slow-moving neighbors. 



B. variable speed model 

Our model system consists of polar particles moving 
in a plane with periodic boundaries. At every time-step, 
each particle updates its position according to the rule 



n{t + At) = nit) + v,{t)h,{t). 



(2) 



To compute the new position we must first define update 
rules for the speed and direction. The latter is given by 



n^{t+l) = 
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(3) 



were Wi is a normalization factor chosen so that fii is a 
unit vector. Noise is introduced by adding the randomly 
oriented vector fji = rj [cos 4>i,sva4>i), with -q the noise 
intensity and a uniformly distributed random variable 
in the interval [— tt, -I-tt]. 

Guided by the experimental results, we formulate a 
minimal variable speed model by considering a simple 
power-law relationship between the local order Xi around 
particle i and the particle speed v^. We write it as 



(4) 



where Mi is the number of individuals in Ui 



If all neighbors within the interaction range f of particle 
i have a similar heading, we will have Xi — ^ ^J^d the 
particle speed will be close to its maximal value vm- By 
contrast, in disordered regions — and particles will 
advance with speed close to zero. Note that for any 7 
an isolated particle will move with maximal speed vmi 
since a particle is always contained within its own neigh- 
borhood, i £ Ui, which implies Xi = 1 in this case. The 
exponent 7 controls the shape of the curve that relates 
local order and speed, as shown in Fig. ^h). For 7 = 0, 
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FIG. 1. (A): Snapshot of experimental system with 300 
golden shiners swimming in a shallow tank. The white circle 
defines the local neighborhood Ui of the focal fish i marked by 
the black dot. The full experiment consists of 3 x 56 minutes 
of video at 30 fps. (B): Experimental relationship between 
the local polarization order x around an individual fish and 
its speed v. The contour lines delineate a 2d surface obtained 
by measuring the probability distribution of v for every given 
value of X- The overlaid curves display the mode of these 
distributions, together with the v = wmX^ relationship used 
in our variable speed rule, for different values of 7. We set 
here vm = 12 in order to make v coincide at x = 1 with the 
mode of the corresponding experimental distribution. 



we recover the fixed speed model. Other examples of the 
variable speed rule in Eq. Q for different values of 7 are 
overlaid on Fig.[TJb). We observe that the variable speed 
rule greatly simplifies the individual dynamics in that it 
replaces the speed distribution expected for a given Xi 
by a single imposed speed value. It captures, however, 
the qualitative dependency of the typical (most common) 
value of Vi on the local polarization order. If we would 



consider the mean of the speed distribution, instead of 
its mode, a similar curve would be obtained but now Vi 
would not approach zero for low Xi values. Either way, 
our minimal variable speed rule Q was not formulated to 
capture these subtleties and provides a reasonable simple 
qualitative approximation of the actual speed dynamics. 

Keeping N and f fixed, the main parameters of our 
model are the mean number density of particles (per in- 
teraction zone area) = Nirf'^ /L"^, the noise intensity 
?7, and the variable speed exponent 7. In our numerical 
study, we will explore different ps regimes by varying the 
system size L and study the order-disorder transition as 
a function of 77 and ps- 



III. ANALYTICAL RESULTS 

Systems of active particles are by construction out of 
equilibrium, since energy is being continuously injected 
at the particle level for self-propulsion. Therefore, free 
energy functional methods cannot be applied. It is, how- 
ever, possible to study the coarse-grained dynamics of 
self-propelled particles in terms of a set of hydrodynamic 
equations, which are derived either using symmetry argu- 
ments |6l |8] or directly from the underlying microscopic 
model [n] [I2I [23j . We choose here the latter approach 
to obtain hydrodynamic equations for the variable speed 
model. Specifically, we will find equations for the coarse- 
grained local density and polarization order parameter 
fields and compare them to the constant speed case |12j . 
The derivation is carried out in the spirit of the Ginzburg- 
Landau theory 24 , that is, in an approximation that as- 
sumes slow modulations in time and space of the local 
order parameters considered. We follow here the same 
approach detailed in |12| . but applying it to a variable 
speed case. 

We begin by defining the coarse-grained local density 
field as 



N 



the polarization order parameter field as the vector 



P(f,i) 



and the traceless symmetric tensor. 



p{r,t) 



(5) 



(6) 



(7) 



which would correspond to the ncmatic order parameter 
field in systems with apolar order [T^]. Here, 5{-) is the 
Dirac delta function, 1 is the identity matrix, and the 
outer product fiihi yields a 2 x 2 matrix with entries 
nfn^, where nf and nf are the a and /3 components 
of unit vector n^, respectively. Using the update rules in 
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Eqs. ([2|, ([3| and Q, together with the analysis in [55], we 
wiU find a stochastic partial difltercntial equation for the 



dynamics of the coarse grained density field p{r,t). We 
start by performing the second order Taylor expansion 



N 



p{r,t + At) - pif,t) = ^ [<5(f- r,(t + At)) - Sir- n{t))] 
1=1 

N ^ N 

W W • ^S{f- r,{t)) + ^Y. W fi^{t)h,{t):VV5{f- n{t)). 



(8) 
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Here, the operator ':' is the double-dot (or colon) prod- ing Vi{t) from the variable speed expressions Q and ([T]), 
uct defined by ab: cd^ J2a J2p a"b^c°'d'^, with indexes dividing by At and taking the fimit A< 0, we find 
a and /3 indicating the vector components. By replac- 
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(9) 



We can write this expression in terms of p, P and Q by 
using definitions ^ and 0. The resuhin g partial 
differential equation describes the dynamics of the den- 
sity field in the Ginsburg-Landau approximation. It is 
given by 



dt 



^7/2 



2 rrP 



(Q+ 2i)p, 



(10) 



where m(r, t) is a coarse-grained field representing the 
number of particles within a circle of radius f and cen- 
tered at r. More precisely, fh{f,t) — ■nf'^^ p{r.j,t), with 
the sum performed over all j satisfying \ fj — 



We are now in a position to compare Eq. ( 10 ) with the 



< f. 



corresponding expression obtained in |12j for particles 
moving with constant speed vc'- 



dp 
di 



-«cV • (Pp) 



(Q + -l)p. (11) 



Note that Eqs. ilOn and (111 have the same form, and 



become identical when we replace 



vc{r,t) 



m(f*, t) 



7/2" 



(12) 



Hence, within the current approximation, we find that 
the particle speed and coarse-grained local number of 
neighbors will be correlated through an inverse power- 
law. This amounts to a relationship between local speed 
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and local density that results from the imposed depen- 
dency of the particle speed on the local polarization or- 
der. We confirm in Section |IV| that this relationship is 
satisfied in numerical simulations, within the disordered 
regime. 

The result in Eq. combined with the variable 

speed rule in Eqs. ([T]) and Q, implies that regions of 
higher density will tend to be less ordered in this regime. 
This is in qualitative agreement with what we have ob- 
served experimentally, but is opposite to the typical rela- 
tionship between mean density and global order found in 
minimal self-propelled particle models TUIIT^. Indeed, 
we will show numerically in Section [Tv| that , even in our 
current variable speed case, the level of polarization or- 
der of the whole system, -0 (defined in Eq. (18) below), 
grows with its mean density ps- Despite its simplicity, 
our model is therefore able to capture a mechanism that 
relates local levels of density and order nontrivially if we 
impose the current particle speed rule. 

By following the same procedure as outlined above (see 
[12 for details), we can also write equations of motion 
for the polarization order parameter. After a long but 
straightforward calculation, we find 



djPp) 
dt 



= F + G + H, 



(13) 



where F is the polynomial term, G the derivative term. 



5 



and H is the noise term. The polynomial term is given by 



F 
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(14) 



The derivative term is 



G 



PV • (Pp) + V(|P|2p) + (p . v)Pp + V • (Qp) + - Vp 



+ 



4mT 



T + V2(Pp)+V(V-P)p . (15) 
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where we have simplified notation by introducing the vec- 
tor f , with T, ^ pViVkiQikP; + P^PfcP^ + PlQ^k) and 
indexes j, fc, and ^ labeling the corresponding tensor com- 
ponents. Finally, the noise term (in the Ito interpreta- 
tion) is found to be 



H = y^Mh. 



(16) 



Here, h{f,t) is a vector field of unit length and random 



while 



orientation, delta correlated in space and time, 
M(r, i) is a 2 X 2 tensor field satisfying = 1. 

We now compare Eq. ( 13 ) to the constant speed case 
derived in (12j. First, we find that the expressions for F 
and H remain unchanged. This implies that the transi- 
tion point, which can be computed in both cases using a 
mean field approximation, will be the same. We confirm 
this result in Section [IV| through numerical simulations. 
By contrast, the derivative term (151 differs from the 
constant speed case, where we have 



G=- 



vc\/m 



PV • (Pp) + V(|P| V) + {P ■ ^)Pp + V • (Qp) 



2^^ 
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r + v2(Pp) + v(v-p)p . (17) 



Again, we find that the constant and variable speed cases 
are equivalent if we replace vc using Eq. (12). 



We conclude that the hydrodynamic equations for the 
density and order parameter fields can be obtained in the 
variable speed case by simply replacing Eq. ( 12 ) into the 
constant speed expressions. 



IV. NUMERICAL STUDY 

In this section, we use numerical simulations to ex- 
plore the range of validity of the analytic description de- 
rived above and study the dynamics of our variable speed 
model beyond this regime. 

We implemented agent-based simulations of = 2000 
particles in a two-dimensional periodic box of side L us- 
ing Eqs. ([2]), ([3]) and All runs presented in this Paper 
were carried out for an interaction range f = 2.0, maxi- 
mum particle speed vm = 0.1 and time-step Ai = 1, so 
that the particle displacement per time-step was never 
greater than 1/20-th of the interaction radius. We stud- 
ied simulations for different levels of noise rj and mean 
density p^, using the variable speed rule in Eq. ^ with 
exponent ranging from 7 = to 7 = 6. The mean den- 
sity was varied by changing the box size L while keeping 
the total number of particles fixed. 

We characterize the collective dynamics resulting from 



simulations using two different global order parameters. 
First, the global polarization order tp, given by 



N 



1=1 



(18) 



This measures the degree of alignment in the system. It 
is equal to 1 when all particles are heading in the same 
direction, and to when they are randomly oriented, re- 
gardless of the particle speeds. Second, the mean particle 
speed u, defined as 



1 ^ 



(19) 



This order parameter is equal to when all particles are 
immobile and to vm when they advance at their maximal 
speed. Note that both quantities are defined here for a 
specific simulation snapshot at time t. We will use below 
these instantaneous values and their averages over time: 
(-0) and (u). 



A. Range of validity of analytical results 

In order to verify the validity of the analytical deriva- 
we confirmed that the relation in 



tions in Section III 
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FIG. 2. Relationship between the particle speed and local 
density that emerges from simulations of our minimal variable 
speed model for 7 = 6. The top and bottom error bars show 
the standard deviations computed using only speeds above or 
below the mean, respectively. The dashed line displays the 
V = vmJ^''^^ relationship deduced analytically in the con- 
tinuous approximation. As the noise is increased and local 
ordered structures vanish, this approximation becomes more 
valid. 



Eq. ( 12 ) holds approximatively for a range of simulation 
parameters. We performed runs using 7 = 6, ps = 2.5 
and various levels of noise rj. We then measured the speed 
Vi of each particle and the number of neighbors A/i within 
its interaction range (including itself). 

The average speed of all particles with a given number 
of neighbors Af is plotted on Fig. [2] In regimes where 
our analytical approach is valid, Eq. ( |12[ ) implies that 
this speed should be approximately equal to vmIM^I"^. 
The figure shows that this relationship is satisfied for high 
levels of disorder, where no ordered structures that could 
violate the approximations in our Ginzburg-Landau ap- 
proach are present. We find that for r\ > 0.7 the in- 
verse power-law proportionality is already verified and 
that above rj 0.85 the exact analytical dependency in 
Eq. ( 12 ) is approximately followed (displayed on the fig- 



ure as a dashed line). The analytic approximation fails, 
however, for noise values < 0.7, where particles are more 
ordered. 



B. The order-disorder transition 

It is common to find in self-propelled particle models, 
such as the Vicsek model [5 , a non-equilibrium phase 
transition that separates the disordered state where par- 
ticles move in random directions from the ordered one 
where they have a common heading. The ordered phase 
can be reached by decreasing the noise level or by in- 
creasing the mean density [9l [26H29] . This transition is 
also present in a variation of the Vicsek model introduced 
by Gregoire and Chate in that is almost identical to 
our current model, but with constant particle speed (i.e.. 
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FIG. 3. Mean global polarization order (tp) (a) and mean 
speed (v) (b) as a function of the noise strength rj for different 
values of the mean density ps- At Tjc « 0.63, the system 
undergoes a discontinuous transition from an ordered to a 
disordered state. This value appears unchanged for the 7 = 
case (constant speed) and for ps > 0.98. As the system 
approaches the disordered state, {v} decreases, vanishing for 
■q > rjc. All curves were computed using = 2000 particles, 
interaction radius f = 2.0, maximum particle speed vm = 0.1, 
and variable speed parameter 7 — 6.0. The different densities 
were obtained by changing the system size, with L — 27, 40, 
80, 120, and 200. 



our 7 = case) . We will examine below the effect of the 
variable speed rule on this transition. 

Figure [3] displays the mean polarization order param- 
eter (f/;) and particle speed (v) as functions of the noise 
ij for different values of ps- We observe that for values of 
Ps >~ 1 there is a sharp discontinuous transition point 
at rjc « 0.6. For lower values of ps, the transition ap- 
pears smoother and, as expected, occurs at lower critical 
noise levels. The decrease in the order parameter (■(/;) is 
accompanied by a reduction of the mean particle speed. 
Note that (v) is already substantially reduced for ry val- 
ues where ('0) is still high. In addition to the curves 
obtained with the variable speed model, we also include 
in Fig. [sf^a) a curve displaying the transition for a con- 
stant speed case (7 = 0) with ps — 8.61. We observe 
that the transition point does not change strongly, as 
predicted by the analytical results presented in Section 
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FIG. 4. Mean global polarization order (tp) (a) and mean 
speed (v) (b) as a function of the mean density ps for differ- 
ent levels of noise ry. All other parameters are the same as in 
Fig. [3] As Ps is decreased, the system undergoes a transition 
from an ordered state with high (tp) and (v) values to a dis- 
ordered state with (ip) « (v) ~ 0. The transition displayed 
by the order parameter (tp) is very similar for the constant 
speed case displayed (7 = 0). For 77 = 0.7 > 77c no transition 
is observed, since the stationary solutions remain disordered 
for all values of the density ps- 



mil We verified that this is also the case for other values 
of 7 between and 6 (data not shown). 

Figure |4] shows the same order-disorder transition as 
Fig. [3l but now as a function of the mean density ps , for 
different values of 77. We find again that the mean particle 
speed decreases as the system loses order. For rj > 0.6, 
the dynamics remain disordered even at high density val- 
ues, which is consistent with the results in Fig.|3] As the 
noise level is reduced, the critical ps decreases until it 
reaches a point where no transition is observed even for 
very low mean denisity values. Finally, we include again 
in Fig. |4|^a) a curve for 7 = to show that the transition 
point also remains unchanged here with respect to the 
constant speed case. 

We conclude from this analysis that the main features 
of the order-disorder transition remain unchanged in the 
variable speed predicted by our analytical calcu- 

lations. However, the critical slowdown of particles that 
is shown above to occur as the system loses order has sig- 



nificant effects in the resulting collective dynamics. We 
will study this phenomenon in more detail below. 



C. Bi-stable solutions 

We now take a closer look at the order-disorder tran- 
sition. Our aim here is not to extrapolate its properties 
to the thermodynamic limit, as this would require larger 
computations and a systematic finite size scaling analy- 
sis. Instead, we focus on understanding how the variable 
speed affects the particle dynamics near the transition. 

Figure [sjj a) displays the global order parameter ip sls a. 
function of time for three different values of 77 close to the 
critical noise. The corresponding histograms are shown 
in Fig. |6] We observe that the dynamics is bistable near 
the transition; the system switches between ordered and 
disordered states. This is in agreement with the first- 
order transition that had been previously reported in the 
constant speed case, both in simulations and in the mean- 
field approximation [5] I26tl29l . We confirm in our simula- 
tions that the 7 = case also displays bistable dynamics 
near the critical point, as shown in Fig. [5jc). We find, 
however, that there is a strong difference between the 
time series obtained in the variable speed and constant 
speed cases. While the former presents very few transi- 
tions between states, even for the long time-series pre- 
sented, the latter shows several intermittent jumps from 
one state to the other [3D]. This difference is explained 
by the lower speed imposed in our model to particles 
surrounded by low local order, which results in a critical 
slowdown of the global dynamics as the order parame- 
ter ip fluctuates to its lower metastable state. Indeed, 
Fig. [sjb) displays the mean particle speed v{t) for the 
same simulation used to generate Fig. [sja) . The corre- 
sponding curves are strongly correlated, with the average 
global dynamics following the local variable speed rule. 
The disordered branch of the bistable solution therefore 
acts similar to an absorbing state; once it is reached the 
particle dynamics slow down so much that it is very hard 
to escape. 



D. Spatial dynamics and phase segregation 

We turn our attention to the spatial dynamics observed 
in the variable speed model. This is where the strongest 
differences with the constant speed case emerge. We find 
that, in the disordered state close to the transition, a 
phase separation occurs where some particles condense 
into almost immobile high-density clusters while the rest 
form a low-density ballistic gas between them. 

Figure [7| presents simulation snapshots of the variable 
speed model with 7 = 6, = 2000, ps = 8.61 and three 
different levels of noise. The length of each line is pro- 
portional to the speed of the corresponding particle. The 
four panels display: an ordered state (a), a disordered 
state (b) , and the two branches of the bistable state dis- 
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FIG. 5. Global polarization order xf}{t) and mean speed v{t) per frame, as a function of time. The variable speed (a and b) 
and constant speed (c) cases are displayed for the same parameters used in Fig.jsj with ps — 8.61 and different noise strengths, 
close to the transition point. Bistable solutions are observed in both cases, but in the variable speed case the transition to the 
disordered branch is accompanied by a critical slowdown of the dynamics, due to the imposed relationship between particle 
speed and local polarization order. The intermittent switching between states thus becomes much less frequent. Note that the 
critical noise value 77c is slightly higher in the variable speed case. 
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FIG. 6. Distribution of the global order parameter ip for the 
variable speed case displayed in Fig.jsja) . The system exhibits 
a discontinuous transition with two coexisting states near the 
critical point. 



cussed in the previous section (c and d). Snapshot (a) 
corresponds to a simulation with low noise i] = 0.1. The 
system is ordered (77 is much lower than its critical value 
r]c ~ 0.625), with particles moving in the same approx- 
imate direction. Since local neighborhoods are mostly 
ordered, x takes values close to 1 in Eq. Q and parti- 
cles advance at speeds close to vm- By contrast, snap- 
shot (b) corresponds to a high noise, r/ = 0.7, disordered 



case. Here, particles are almost immobile. Due to the 
high mean density, even if the noise spreads out particles 
almost homogeneously, each one has several neighbors 
within its interaction range. Given the low local order, 
X ~ in Eq. Q and particle speeds are close to zero. 

Snapshots (c) and (d) are particularly interesting as 
they display previously unobserved dynamics, not possi- 
ble in a constant speed model, that occur when the noise 
is close to its critical value. As shown in Figs, [s] and 
[6) we have in this case bistable global dynamics. The 
two metastable solution branches correspond to an or- 
dered state with high and to a disordered state with 
low ■0. Both states organize into a high-density cluster 
surrounded by a low-density particle gas. In panel (c), 
the cluster moves as indicated by the red arrow, while in 
panel (d) it remains almost static. Larger simulations can 
display several clusters with similar behavior. Particles 
surrounding the clusters move rapidly. 

While the presence of moving clusters is common in 
Vicsek-like models [301 IS], static clusters such as the 
one on Fig. have not been previously observed. The 
mechanism that leads to their nucleation therefore de- 
serves more a detailed analysis. As the system fluctuates 
between an ordered and a disordered state (being close 
to the transition point), it reaches situations where par- 
ticles are typically not locally aligned. This leads to the 
formation of groups of slow moving particles that grow 
as other particles reach them, because incoming parti- 
cles suddenly confront high density disordered regions, 
thereby losing orientation and slowing down drastically. 
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FIG. 7. Simulation snapshots for — 2000 particles in four different steady states with the same mean density pa — 8.61, 
maximum particle speed vm ~ 0.1, and variable speed exponent 7 = 6. Each particle is represented by a short line with 
length proportional to its speed, (a) Low noise {r] = 0.1) ordered state. We observe clusters of moving particles (with velocity 
indicated by the arrow) equivalent to those found in models with constant speed, (b) High noise {i] — 0.7) disordered state. 
Particles are almost immobile and homogeneously distributed in space. Due to the high ps value considered here, each particle 
interacts with several neighbors heading in different directions, resulting in low particle speed, (c and d) Intermediate noise 
states (77 = 0.6253) close to the transition. The system here can be either in the ordered solution branch (c), displaying a moving 
cluster (with velocity indicated by the arrow), or in the disordered branch (d), which exhibits a static cluster as described in 
the main text. 



Once clusters are formed, they can only lose particles at 
their frontier. There, some particles will spontaneously 
align due to noise fluctuations and manage to escape 
if they head away from the cluster, since they will feel 
less and less the influence from the disordered region. 
Eventually a single particle, or a small group, will be 
far enough removed to feel only its own influence, form- 
ing a low density ballistic gas between clusters. Particles 
in this gas are isolated and their x is therefore close to 
1. They move at speeds close to vm until they reach 
another cluster and condense again. These clusters will 
thus grow until their absorption rate is balanced by their 
evaporation rate. In the bulk of the clusters particles 
advance very slowly, which hinders their ability to re- 
order. Indeed, it has been shown that, when noise is 
present, Vicsek-like models can only reach an ordered 
state if particles are able to move with respect to each 
other, which allows them to switch neighbors and estab- 
lish effective long-range interactions over time [221 . If 
particles always interact with the same neighbors, the 
system becomes equivalent to an XY-model, for which 
the Mermin- Wagner theorem shows that no long-range 
order can exist in two dimensions at nonzero noise levels 
|32j . Hence, the lack of relative particle motion in the 
bulk of a static cluster helps stabilize it in a disordered, 
immobile state. 

In order to study how the particles distribute between 
those that are part of a static cluster and those that are 
not, we plot on Fig. [8] the individual particle speed dis- 
tribution for various values of noise and density. Panel 
(a) shows that, as the noise level is incresased, first the 
typical particle speed is smoothly reduced and then a 
zero-speed peak emerges, corresponding to static clus- 
ters. While in the ordered phase, groups of particles 
therefore move slower and slower until the critical noise 
level is reached and static clusters start nucleating. In 



panel (b) the noise level is fixed at a high value {rj = 0.7), 
so the system is always in the disordered state. As we 
vary the density ps, the particle speed distribution re- 
mains bimodal. Particles can either be within a static 
cluster and stop advancing or be part of the inter-cluster 
gas and move at maximum speed vm while remaining iso- 
lated. As the mean density is increased, more particles 
are trapped in static clusters and thus have zero speed. 

FinalUy, we present in Fig. [9] the individual particle- 
speed distribution as a function of the variable speed 
exponent 7. We consider two parameter regimes, one 
with noise intensity just below the critical noise value 
[Fig. [9]ja), 77 = 0.6] and one just above it [Fig. ^h), 
77 = 0.7]. In the first situation, the typical particle speed 
gradually decreases as 7 is increased, from v ss vm = 0.1 
when 7 = 0, to u sa 0.04 for 7 = 6. In the disordered 
state, by contrast, the typical particle speed is reduced 
much more abruptly for any 7 > 0, quickly giving rise to 
the nucleation of static clusters. 



V. DISCUSSION 

The analysis above demonstrates the similarities and 
differences between standard, minimal constant speed 
models of collective motion [51 [TO] and our variable speed 
version. We show that the variable speed case displays 
an order-disorder transition analogous to that observed 
for constant speed. In both cases, as the noise level is de- 
creased or the mean density increased, the system goes 
from a disordered state to an ordered one where they 
align to a common heading. The dynamics associated to 
the transition, however, are very different. 

An interesting finding is that our variable speed rule 
induces an inverse power-law relationship between the 
speed of a particle, or the level of order in its surround- 
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FIG. 8. (a) Distribution of individual particle speeds for the 
same parameters as in Fig. |3j with ps = 8.61 and different 
levels of noise intensity. As the noise level is increased, the 
typical particle speed decreases until static clusters start nu- 
cleating for rj > rjc ~ 0.625. These clusters absorb most par- 
ticles, while a few isolated particles continue to move between 
them at maximal speed, as shown by the peaks at u = and 
V ~ vm = 0.1. (b) Distribution of individual particle speeds 
for the same parameters as in Fig. [4] with ri — 0.7 (disordered 
state) and different values of the mean density ps. Here, static 
clusters are always present, their size increasing with ps- Par- 
ticles are again either frozen within these clusters or moving 
alone between them at maximal speed. 



ings, and the local density. Note that this is opposite 
to the typical correlation between global order and mean 
density in minimal models [5] . Indeed, high mean density 
makes particles interact in average with more neighbors, 
which increases the level of order, as it has been shown 
numerically and analytically (in the mean field approxi- 
mation) 28, 29 . 

While the model presented in this paper is far from re- 
alistic, the inverse relationship between local quantities 
specified above appears to be robust enough to persist in 
less idealized systems. In the golden shiner data we find 
not only the correlation between individual speed and 
local polarization order that inspired our variable speed 
rule (see Fig. [T]) , but also an inverse relationship between 
local order, or particle speed, and local density analo- 



FIG. 9. Distribution of individual particle speeds for the 
same parameters as in Fig. [3] with noise intensity just below 
{rj — 0.6) or above (rj = 0.7) the critical noise value, and 
different values of the variable speed exponent 7. When in 
the ordered state (a), the typical particle speed gradually de- 
creases as 7 is increased, from a fixed v = vm = 0.1 value for 
7 = 0, to v « 0.04 for 7 = 6. In the disordered state (b), 
by contrast, the typical particle speed is reduced much more 
abruptly for any 7 > 0, where the static clusters described in 
the main text are nucleated. 



gous to that deduced from our model (data not shown). 
We emphasize that we cannot deduce a causal relation 
between these effects from the experimental data. By es- 
tablishing connections between them in a specific model, 
however, we have provided a mechanism through which 
they could be robustly related in generic situations. 

Another observation from our study is the nucleation 
of static clusters. In simulation videos the dynamical pro- 
cess leading to these clusters looks similar to jamming in 
granular materials (33J. Despite the differences between 
these processes, we can draw some analogies that go be- 
yond their visual appearance. In physical systems, the 
opposing forces on jammed particles add to zero, stop- 
ping their flow. Likewise, here opposing headings (i.e. 
conflicting information) add up to zero, producing no lo- 
cal order and therefore a vanishing particle speed. We 
also note that both jammed and static regions will grow 
by recruiting moving particles that reach them and are 
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stopped by opposing interactions. Finally, an additional 
similarity is given by the dynamics in the bulk and edges 
of the clusters. In both cases edge particles can escape if 
the sum of all surrounding interactions points away from 
the cluster, while bulk particles can only move if there is 
a pathway of particles with a non-zero sum of interactions 
along it that percolates through the group [55] . 

Despite the similarities outlined above, the mecha- 
nisms behind both processes are very different. Jam- 
ming is produced by contact forces or repulsive poten- 
tials, while interactions in the variable speed model are 
based on heading directions. Another important differ- 
ence is that in our model static clusters cannot form be- 
low a certain noise level, while jamming is always in- 
creased at lower noise levels. 

We end this section by pointing out that our variable 
speed model considers only one possible way of relating 
the speed of an individual particle to its local environ- 
ment. While the model was inspired by experimental 
data, we do not claim any specific causal origin to this 
correlation. It could result from the interplay of a num- 
ber of biological interactions that are not considered in 
this minimal model. It could also be for example, that 
alignment is enhanced in faster moving particles or that 
the speeding rule depends in fact directly on the local 
density. 

VI. CONCLUSION 

We have studied the dynamics of a minimal model of 
collective motion where the particles move with variable 



speed. We found that, despite the simplicity of the al- 
gorithm and its similarity with standard constant speed 
models, our system displays not only the usual ordering 
transition but also new generic dynamical phenomena 
that could be present in more realistic models or in ex- 
perimental systems. These are: (1) an inverse correlation 
between individual particle speed or local polarization or- 
der and local density, produced by the imposed relation- 
ship between particle speed and local order, and (2) the 
nucleation in certain regimes of static clusters where in- 
dividuals remain immobile, only turning in place without 
achieving order. 

We conclude that including variable speed dynamics in 
standard self-propelled particle models produces a range 
of new phenomena that can be relevant for experimental 
systems. Given the amount of work dedicated up to now 
to models with constant speed, we hope that our work 
will serve as a motivation to explore the generic and spe- 
cific consequences of considering variable speed rules in 
this class of systems. 
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